Exact final state integrals for strong field QED 
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This paper introduces the exact, analytic integration of all final state variables for the process of 
nonlinear Compton scattering in an intense plane wave laser pulse, improving upon a previously slow 
and challenging numerical approach. Computationally simple and insightful formulae are derived for 
the total scattering probability and mean energy-momentum of the emitted radiation. The general 
form of the effective mass appears explicitely. We consider several limiting cases, and present a 
quantum correction to Larmor's formula. 
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Introduction: The recent and predicted progress in 
laser technology leading to very high peak intensities jus- 
tify the need for a better understanding of so-called non- 
linear QED, describing phenomena occurring in fields so 
strong that their effects cannot be treated perturbatively. 

Unfortunately, the complexity of the processes inside 
these ultra-intense laser beams has meant that several 
simplifications have had to be used to make practical 
computations feasible. The laser beam is usually sup- 
posed to be in a coherent state, which can be well ap- 
proximated by a classical field. Due to the relatively 
small frequencies, unless massive particles of very high 
Lorentz factor are involved, quantum effects are small, 
so a fully classical description may often be justified. For 
instance, in discussing the scattering of an electron in a 
laser beam, one may consider Thompson scattering in- 
stead of its quantum counterpart, nonlinear Compton 
scattering (NLCS) fl. The classical approximation al- 
lows for a realistic description of the laser field and the 
inclusion of radiation reaction (RR) 2], but is unable to 
describe important quantum effects, such as nonpertur- 
bative pair creation from vacuum Q , the trident process 
or vacuum birefringence 

A treatment of these processes in the framework of 
nonlinear QED, even in a semiclassical approach, has 
not yet been performed without further approximations, 
such as replacing the laser field by an idealized plane 
wave, thus allowing for analytical (Volkov) solutions to 
the Dirac equation. This disregards the strong spatial fo- 
cusing of the beam needed to attain high intensities. In 
addition, for a long time, results were restricted to only 
infinite, monochromatic plane waves, or even, giving up 
periodicity, to a crossed field model [6, 7]. Only recently, 
the more realistic short pulse plane waves came into use 



In [i, 11, l3|, the photonic and electronic distributions 
resulting from NLCS were described in detail for some 
model pulses. The total probability and emitted parti- 
cle's energy-momentum can be obtained by integrating 
these distributions, but the task is numerically challeng- 
ing. By a change in integration order and a different 
regularization, allowing for all final state integrals to be 
performed analytically, we obtain equivalent, but more 
revealing and much easier to compute formulae. For sim- 



plicity, we only present unpolarized results. The method 
we develop is quite general and shall be applied to other 
strong field QED processes in a future paper. 

Preliminaries: For our starting point, all notations 



and conventions, we refer the reader to 13|. In short, 
p — mv and k' are the initial electron and final photon 
four-momenta. We opt for natural units, so c = = 1. 
Define k = ojn, where — and oj is some characteristic 
frequency of the wave. Let (j) = k ■ x — ujx~ be an invari- 
ant lightfront coordinate, used to describe the plane wave 
pulse by the four-potential A — ^aof (0). By transver- 
sality, k ■ A — 0. We choose n = (1,63), /o = /a = 0, 
and use lightfront notations such as = p^ ±p^ , p-*- — 
p — p'^n. The final results for probability/momentum 
will prove to be manifestly Lorentz invariant /covariant. 
While working with / ((/>) , only the gauge changes keep- 
ing A </)-only dependent are allowed. But the end results 
can be expressed in terms of the classical velocity, so they 
are gauge invariant. For a long pulse, one may choose lo 
the carrier frequency and the peak value of the envelope 
of /' ((f)) equal to one. To compare very short pulses one 
may prefer to fix w and ao so that the pulse's range is of 
order 27r and — J d<j>f''^ {(f)) ~ 1. Whatever our choice, oq 
should offer a reliable description of the peak intensity. 
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Classical motion: Let tt = mu be the kinetic momen- 
tum of a classical electron moving in this plane wave field, 
where 



u{(f) = v- aaf {(f)) + 



2aof {^)-v- alP {^) 
2k ■ V 



k. 



(1) 



If we set / (—00) = 0, then v = u (—00) is indeed the 
velocity of the particle before meeting the wave. As op- 
posed to unipolar pulses that permanently accelerate the 
particle [l^, for the usual whole-cycle pulses, / (00) = 
and u (00) = u (— cxd), as long as RR is neglected. 

Scattering probability: In the Furry picture, the total 
NLCS probability for unpolarized initial electrons is: 
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where 9 ~ 



5(0 = 



and we denoted the 



moving average of a function F by (F) 
For instance, (/') = [/ ((/)') - / (0)] /O. The formal ex- 
pression (21) was obtained from formula CI in using 
the unregularized form of — /^e~**c?0, with € 
(l, /i, /2, f^), lightfront coordinates and a 4-dimensional 
integration. See also All B^, in the generic case, 

or at least Bq, need to be regulated, damping the oscil- 
lations of the integrand with a convergence factor such 
as e"^*^ , e \ 0, that can be discarded after a partial 
integration restricts B^, to the length of the pulse 3, [13] ■ 
At first glance, in writing ([2|) we have added to the 
numerical complexity, constructing the inner double in- 
tegral from the simple integrals B,^. But, in fact, by 
a change of quadrature order, the analytical integration 
over k'-*-, and k'~ leaves us with only two easy integrals, 
instead of the initial four. In addition, we get rid of the 
rapid oscillations encountered when computing B^. Ex- 
pressing k' ■ (tt) in lightfront coordinates, we notice the 
integral over k'^ is Gaussian, if regulated by replacing in 
the exponent the factor 9 hy 9 + ie, then taking e \ 0. 
The previous damping factor is now superfluous. Intro- 
ducing the invariant number bn — the result is 
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where 



^i = l 



> 1. 



(3) 



(4) 



We recognize the effective mass M = d (tt) , first intro- 
duced by Kibble, that appears in the Volkov propagator 
14 [ and the Wigner function 15, 16j : 



pi = {uY 



(5) 



A Lorentz and gauge invariant, M depends on the aver- 
aging interval. For a whole-cycle finite pulse, the mass 
shift M — m vanishes when 9 ^ oo. We now use the 
relation {9 + ie)~^ = p.v.9~^ — iirS {9) and the fact that 
H is an even function of 9, so the result is indeed real. 
Changing variables from 0, (j>' to 9 and a = "^^"^ , and 
noticing that tt = Jjj da; sine a; , we get 



P = - 
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A new analytical integration leads to: 
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where, in terms of trigonometric integrals, 

Ji (x) = -xA' (x) , J2{x)^l[2 + x~ x^A (x)] , 

o 

A (x) = sin X ci (x) — cos x si (x) , 
A' {x) — cos X ci (x) -I- sin x si (x) . 

Notice that the Lorentz invariant ([7]) depends on 6o, and 
4-products of the values of the function aof, but not 
on aof ■ p. Interestingly, we can rewrite the probabil- 
ity in terms of the classical velocity ([l} as function of 
the proper time, eliminating all reference to the driving 
field and proving gauge invariance. Could the result be 
generalized to an arbitrary, not necessarily plane, wave? 
It is hard to answer, because of the many very different 
trajectories allowed inside a general field. The attrac- 
tiveness of a plane wave derives from the simplicity of 
the law of motion it entails. The electron's motion al- 
ways looks the same, regardless of the initial position. A 
quantum computation for a general field would require 
the use of a wavepacket with some initial average posi- 
tion and momentum that only in the limit relates to a 
particular classical motion. 

Radiated energy-momentum: The same procedure 
can be applied to compute the expectation value of the 
emitted photon's momentum. 
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Following the same method as for the probability, we are 
left with the manifestly covariant double integral 
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The new functions used are: 

Jg (x) = ^A (x) - xA' (x) - f 

J,{x)^l,[{x'-(Sx^)A{x)]' -^ + \ 

The standard deviation of k"^ or any of its higher mo- 
ments can similarly be computed. Even though derived 
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for a pulse, our formulae can provide the emission prob- 
ability and radiated energy-momentum per cycle in the 
case of an infinite, periodic plane wave, by truncating it 
to only a finite number of cycles N, and considering the 
quickly reached N ^ oo limit of P/N and {k'"} /N. The 
results are given by (0 and © with the a integral re- 
stricted to one period. A simple example is a circularly 
polarized monochromatic wave f (0) — (cos 0, sin0, 0), 
for which /i = 1 + Oq (l — sine is independent of a, so 
J da ^ 2tt. 

Of particular interest are the limiting cases of small gq 
or bo. With today's technology, a; is around the order of 
1 — 10* eV and only for the optical range a large oq is 
possible. If, for instance, a; = 3 eV, an electron energy 
of at least 20 GeV is needed for bo to reach unity. 

Perturbative limit The results can be expressed as 
integrals containing the field's Fourier transform and 
the Klein-Nishima probability of linear Compton [l3| . 
To emphasize the role of 6o, we prefer to work with 
a quadratic function related to the potential's auto- 
correlation. 



one obtains the approximation 
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The weak field limits (ao < 1) of (0 and (O are: 
2a 



P. 



dev 



JO 

dd ' 



.„y^(^). (12) 
' ""f (0), (13) 



where 7^, J-, and Q are the universal, positive functions, 
V{x)^x' lL^J,{y)dy + xJ2{x), 

J X 

T{x)^x^ j y:^J,{y)dy + xJ,{x), 
g{x)^x' ^J^[y)dy + J^{x). 



Since V {x) is increasing for x > 0.fT2l decreases with bo, 
having the classical upper bound 



9 r'^ 
Pp. ci^^ d0£ (9) . 
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Classical limit: This is expected when &o ^ 1, since 
bo is proportional to h, a fact obscured by the use of 
natural units. That is, the laser photon energies are much 
smaller than m in the electron's rest frame. Now the 
major contribution to the integral ^ comes from small 
6. By Taylor expansions, such as 
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where 
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(15) 



1 /"^ 

C{x) = - [6 J4 {v + xy^) ~Jz{v + xy"")] ydy. 
^ Jo 

(16) 

The above function is decreasing for x > 0, quickly de- 
parting from C [0] = 1. Already C [0.01] ~ 0.58. No clas- 
sical equivalent is found for the second term in the inte- 
grand of ([9|), negligible in this limit. An equivalent form 
of (|15p can be written by expressing everything in terms 
of a classical particle's velocity as function of proper time. 
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When not only is 69 ^ Ij but also ^qAq ^ 1, equation 
([T7]) further reduces to: 



le2 /■ 
/ dru" (r) ii (r) 



(18) 



The same result arises directly from classical electrody- 
namics, if one neglects radiation reaction. For z/ = 0, 
([TS]) is just the time integrated Larmor's formula. Notice 
that a small 60 doesn't necessarily make classical Thomp- 
son scattering a good model. For strong fields, much less 
power is radiated than Larmor predicts. In CED, the 
spectral and angular distribution is also a double inte- 
gral, but an interesting cancelling of interference terms 
leaves a single one after the frequency is integrated away. 
The transfer of energy and momentum to the field is well 
defined at each moment in time, there are no quantum 
uncertainties. 

Interestingly, this form of decoherence is also shown 
by the better approximation (I17p . that looks mislead- 
ingly classical, though the argument of the correction C 
is in fact proportional to . Heuristically, bo is a sort 
of coherence lighfront time 9 scale. When small enough, 
it allows for a time-incoherent model of emission, but at 
high intensities the mass shift implied by (ITil) cannot be 
neglected in this coherence interval, hence the aforemen- 
tioned correction. This happens when the peak electric 
field in the electron rest frame approaches the Schwinger 
critical value m? je. Formula ()17|) could provide a general 
improvement to Larmor's, valid in an arbitrary driving 
field, whose frequencies in the rest frame of the electron 
are similarly low compared to m, so the emission can be 
viewed as incoherent in time, hence the product of the 
classical motion of a charged particle. As for the total 
scattering probability, by an asymptotic expansion of the 
functions Ji [x) one gets the limit 
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In this case no decoherence is observed, as photon emis- 
sion probabihty is not a classical concept. 
Numerical example: 
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the shape of the electric field strength 

We present results for a linearly polarized, one cycle, 
pulse, characterized by a Gaussian potential: 



(20) 



The NLCS probability is shown in Fig[TJ The quadratic 
increase from the perturbative region, shown in detail in 
the upper left corner, quickly slows down as oq grows past 
unity. The higher the parameter bo, the lower P is. In 
general, ([7]) boundlessly grows with the length/intensity 
of the pulse. Even for one as short as (|20|) and the ex- 
perimentally attainable oq = 100, the result can easily 
In 



surpass unity. 



this possibility was noticed, inter- 
preted as a sign that multiphoton emission cannot then 
be neglected, and a re-normalization of the whole series of 
n-photon NLCS probabilities was suggested. Moreover, 
for a unipolar pulse, ([7]) shows the logarithmic IR diver- 
gence typical of Bremsstrahlung 13|, llSl j . These problems 
can be dealt with by including a one loop self-energy di- 
agram, that adds nothing to ©, but does contribute to 
the expectation value of the final electron's momentum, 
even in the whole-cycle case [l^ . For the general theory 
of the cancellation between real and virtual photon IR 
divergences, see [20| . 

To plot formula ([9]), we need more than just the in- 
variants oo, b(). Let's assume uj = I eV. It remains to 
know the incidence, that we choose head-on. While ex- 
perimentally difficult, this gives the largest bo for a given 
pulse and electron beam. In Figj^] a comparison is drawn 
between the expectation value of the radiated energy and 
its two incoherent approximations. Both (|18p and PT)) 
overestimate but the latter is a much closer match. 

Conclusions: We have found a way to analytically in- 
tegrate all final state variables out of the NLCS proba- 
bility and expectation values. This not only saves huge 
computational expense, but also reveals new insights into 
the structure of scattering processes in strong fields. We 
shed light on the role of the effective mass and the emis- 
sion's coherence in time. Simple results were found for 
the monochromatic, perturbative and classical limits. We 
derived a strong field correction to Larmor's formula, 
arising from the mass shift. A computation using a real- 
istic model showed the usefulness of this approximation. 
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FIG. 1: The NLCS probability ^ is plotted as function of 
ao, for bo — 10^ , p = —4, 1. The dashed line is the classical 
limit (fT9l). 




FIG. 2: The expectation value of as a multiple of the 

electron rest energy, plotted as function of ao. The dotted 
line is the Thompson limit (|18p . while the dashed one is the 
better approximation (|17|l . 



but also found large values, even surpassing unity, for 
the probability. This suggests multiple scatterings and 
radiative corrections need to be considered. In a future 
paper, our method will be applied to these, as well as to 
other strong field processes. A thorough numerical ex- 
ploration, comparison and interpretation of our results 
for various pulses and setups is in preparation. 
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